Variace konstant
Metoda variace konstant nebo variace parametrů je v matematice obecná metoda řešení nehomogenních lineárních obyčejných diferenciálních rovnic.
U rovnic prvního řádu je obvykle snazší hledat řešení metodou integračních faktorů nebo neurčitých koeficientů. Tyto metody však využívají heuristiky, které zahrnují hádání, a nefungují pro všechny nehomogenní lineární diferenciální rovnice.
Metodu variace konstant lze použít i pro lineární parciální diferenciální rovnice, konkrétně pro nehomogenní problémy u rovnic s lineárním rozvojem, jako je rovnice vedení tepla, vlnová rovnice a rovnice kmitající desky. V těchto případech je metoda známější pod názvem Duhamelův princip, podle Jean-Marie Duhamela, který ji poprvé použil na řešení nehomogenní rovnice vedení tepla. Někdy je přímo variace parametrů nazývána Duhamelův princip a naopak.
Historie
[editovat | editovat zdroj]Metodu variace konstant poprvé použil švýcarský matematik Leonhard Euler (1707–1783) a rozpracoval italsko-francouzský matematik Joseph Louis Lagrange (1736–1813)[1]. Předchůdce metody variace dráhových parametrů nebeských těles se objevil v Eulerově práci v roce 1748, kdy studoval vzájemné poruchy pohybu Jupitera a Saturna[2]. Při studiu pohybu Země v roce 1749 získal Euler diferenciální rovnice pro dráhové parametry[3] a v roce 1753 aplikoval tuto metodu na své studium pohybu Měsíce[4]. Lagrange první použil tuto metodu v roce 1766[5]. V letech 1778–1783 Lagrange metodu dále rozvinul jak v řadě monografií o poruchách pohybu planet[6], tak v další řadě prací o určování drah komet ze tří pozorování[7]. (Je třeba zmínit, že Euler a Lagrange aplikovali tuto metodu na nelineární diferenciální rovnice, a že místo variací koeficientů v lineární kombinaci řešení homogenní rovnice aplikovali variace na konstanty v rovnicích nerušených pohybů nebeských těles[8]). Během let 1808–1810 dal Lagrange metodě variace parametrů její konečnou podobu v řadě prací[9]. Hlavním výsledkem jeho studia byl systém planetárních rovnic v Lagrangeově tvaru, který popisuje rozvoj Keplerovských dráhových parametrů pohybu tělesa rušeného dalšími tělesy.
Ve svém popisu rozvojů drah Lagrange vycházel ze zjednodušeného problému dvou těles pro získání nerušeného řešení a předpokládal, že všechny poruchy pocházejí z gravitačního působení ostatních těles na obíhající těleso. Následkem toho jeho metoda předpokládá, že poruchy závisí pouze na poloze obíhajícího tělesa, ne na jeho rychlosti. Od 20. století se v nebeské mechanice uvažují interakce, které závisí nejen na polohách ale i rychlostech (relativistické korekce, odpor atmosféry, inerciální síly). Proto byla metoda variace parametrů používaná Lagrangem rozšířena na situace se silami závislými na rychlosti[10].
Popis metody
[editovat | editovat zdroj]Je-li dána obyčejná nehomogenní lineární diferenciální rovnice řádu n
nechť je fundamentální systém řešení odpovídající homogenní rovnici
Pak určité řešení nehomogenní rovnice je
kde jsou derivovatelné funkce, tj. předpokládá se, že vyhovují podmínce
začínáme s (iii), opakovaná diferenciace kombinovaná s opakovaným použitím (iv) dává
Poslední diferenciace dává
Substitucí (iii) do (i) a použitím (v) a (vi) dostáváme, že
Lineární systém (iv a vii) n-té rovnice pak lze vyřešit pomocí Cramerova pravidla, což dává
kde je Wronskián fundamentálního systému a je Wronskián fundamentálního systému s i-tým sloupcem nahrazeným
Určité řešení nehomogenní rovnice lze pak psát jako
Příklady
[editovat | editovat zdroj]Zvláštní rovnice druhého řádu
[editovat | editovat zdroj]Máme řešit rovnici
nejdříve najdeme obecné řešení diferenciální rovnice, to jest řešení homogenní diferenciální rovnice
určíme kořeny charakteristické rovnice
protože řešením je vícenásobný kořen, musíme zavést faktor x pro jedno řešení pro zajištění lineární nezávislosti.
Tak získáme u1 = e−2x a u2 = xe−2x. Wronskián těchto dvou funkcí je
Protože Wronskián je nenulový, funkce jsou lineárně nezávislé, takže jsme získali přímo obecné řešení homogenní diferenciální rovnice (nejen jeho podmnožinu).
Hledáme funkce A(x) a B(x), tak že A(x)u1 + B(x)u2 je obecné řešení nehomogenní rovnice. Potřebujeme pouze vypočítat integrály
to jest,
kde a jsou integrační konstanty.
Obecné rovnice druhého řádu
[editovat | editovat zdroj]Máme diferenciální rovnici tvaru
a definujeme lineární operátor
kde D reprezentuje diferenciální operátor. Proto máme jak řešit rovnice pro , kde a jsou známé.
Musíme vyřešit první odpovídajícím homogenní rovnici:
libovolnou technikou, kterou si vybereme. Jakmile budeme mít dvě lineárně nezávislá řešení této homogenní diferenciální rovnice (protože se jedná o rovnici druhého řádu) — označíme je u1 a u2 — můžeme využít metodu variace konstant.
Od tohoto okamžiku tedy hledáme obecné řešení diferenciální rovnice , o němž předpokládáme, že je tvaru
kde a jsou neznámé funkce a a jsou řešení na homogenní rovnice. Všimněte si, že jestliže a jsou konstantní, pak . Požadujeme, aby A=A(x) a B=B(x) byla tvaru
Nyní,
a díky platnosti výše uvedené podmínky dostáváme
Dalším derivováním (mezikroky nejsou uvedeny)
Nyní můžeme zapsat použití L na uG jako
Protože u1 a u2 jsou řešení, pak
dostáváme systém rovnic
roznásobením dostaneme
Takže výše uvedený systém určuje právě podmínku
Hledáme A(x) a B(x) z této podmínky, tak, daný
můžeme řešit pro (A′(x), B′(x))T, tak
kde W označuje Wronskián u1 a u2. (Z předpokladu, že u1 a u2 jsou lineárně nezávislé, víme, že W je nenulový.)
Tedy
Zatímco homogenní rovnice lze řešit relativně snadno, tato metoda umožňuje výpočet koeficientů obecného řešení nehomogenní rovnice a tedy umožňuje zjistit úplné obecné řešení nehomogenní rovnice.
Všimněte si, že a jsou vesměs určeny až na aditivní konstantu (integrační konstanta); očekávali bychom dvě integrační konstanty, protože původní rovnice byla druhého řádu. Přidáním konstanty k nebo se hodnota nezmění, protože je lineární.
Reference
[editovat | editovat zdroj]V tomto článku byl použit překlad textu z článku Variation of parameters na anglické Wikipedii.
- ↑ Viz
- Forest Ray Moulton, An Introduction to Celestial Mechanics, 2nd ed. (poprvé publikováno vydavatelstvím Macmillan Company v roce 1914; reprint v roce 1970 vydalo nakladatelství Dover Publications, Inc., Mineola, New York), stránka 431.
- Edgar Odell Lovett (1899) "Theory of perturbations a Lie's theory of contact transformations," Quarterly Journal of Pure a Applied Mathematics, vol. 30, stránky 47-149; viz zvláště stránky 48-61.
- ↑ Euler, L. (1748) "Recherches sur la question des inégalités du mouvement de Saturne et de Jupiter, sujet proposé pour le prix de l'année 1748, par l’Académie Royale des Sciences de Paris" [Zkoumání poruch pohybu Saturna a Jupitera; toto téma bylo navrženo Královskou akademií věd v Paříži pro cenu roku 1748] (Paříž, Francie: G. Martin, J.B. Coignard, & H.L. Guerin, 1749).
- ↑ Euler, L. (1749) "Recherches sur la précession des équinoxes, et sur la nutation de l’axe de la terre," Histoire [nebo Mémoires ] de l'Académie Royale des Sciences et Belles-lettres (Berlin), stránky 289-325 [publikováno v roce 1751].
- ↑ Euler, L. (1753) Theoria motus lunae: exhibens omnes ejus inaequalitates … [Teorie pohybu Měsíce: demonstrace všech jeho nepravidelností… ] (Saint Petersburg, Russia: Academia Imperialis Scientiarum Petropolitanae [Imperial Academy of Science (St. Petersburg)], 1753).
- ↑ Lagrange, J.-L. (1766) “Solution de différens problèmes du calcul integral,” Mélanges de philosophie et de mathématique de la Société royale de Turin, sv. 3, stránky 179-380.
- ↑ Viz
- Lagrange, J.-L. (1781) "Théorie des variations séculaires des élémens des Planetes. Premiere partie, … ," Nouveaux Mémoires de l'Académie Royale des Sciences et Belles-lettres (Berlin), stránky 199-276.
- Lagrange, J.-L. (1782) "Théorie des variations séculaires des élémens des Planetes. Seconde partie, … ," Nouveaux Mémoires de l'Académie Royale des Sciences et Belles-lettres (Berlin), stránky 169-292.
- Lagrange, J.-L. (1783) "Théorie des variace périodiques des mouvemens des Planetes. Premiere partie, … ," Nouveaux Mémoires de l'Académie Royale des Sciences et Belles-lettres (Berlin), stránky 161-190.
- ↑ Viz
- Lagrange, J.-L. (1778) "Sur le probleme de la détermination des orbites des cometes d'après trois observations, premier mémoire" (K problému určování drah komet ze tří pozorování, první monografie), Nouveaux Mémoires de l'Académie Royale des Sciences et Belles-lettres (Berlin), stránky 111-123 [publikováno v roce 1780].
- Lagrange, J.-L. (1778) "Sur le probleme de la détermination des orbites des cometes d'après trois observations, druhý mémoire", Nouveaux Mémoires de l'Académie Royale des Sciences et Belles-lettres (Berlin), stránky 124-161 [publikováno v roce 1780].
- Lagrange, J.-L. (1783) "Sur le probleme de la détermination des orbites des cometes d'après trois observations. Troisième mémoire, dans lequel on donne une řešení directe et générale du problème.", Nouveaux Mémoires de l'Académie Royale des Sciences et Belles-lettres (Berlin), stránky 296-332 [publikováno v roce 1785].
- ↑ Michael Efroimsky (2002) "Implicit gauge symmetry emerging in N-body problem of celestial mechanics," stránka 3.
- ↑ Viz
- Lagrange, J.-L. (1808) “Sur la théorie des variations des éléments des planètes et en particulier des variace des grands axes de leurs orbites,” Mémoires de la première Classe de l’Institut de France. Reprinted in: Joseph-Louis Lagrange s Joseph-Alfred Serret, ed., Oeuvres de Lagrange (Paris, France: Gauthier-Villars, 1873), vol. 6, stránky 713-768.
- Lagrange, J.-L. (1809) “Sur la théorie générale de la variation des constantes arbitraires dans tous les problèmes de la méchanique,” Mémoires de la première Classe de l’Institut de France. Reprinted in: Joseph-Louis Lagrange s Joseph-Alfred Serret, ed., Oeuvres de Lagrange (Paris, France: Gauthier-Villars, 1873), vol. 6, stránky 771-805.
- Lagrange, J.-L. (1810) “Second mémoire sur la théorie générale de la variation des constantes arbitraires dans tous les problèmes de la méchanique, … ,” Mémoires de la première Classe de l’Institut de France. Reprinted in: Joseph-Louis Lagrange s Joseph-Alfred Serret, ed., Oeuvres de Lagrange (Paris, France: Gauthier-Villars, 1873), vol. 6, stránky 809-816.
- ↑ Viz
- Michael Efroimsky (2005) "Gauge Freedom in Orbital Mechanics." ANYAS, Vol. 1065, pp. 346–374 (2005)
- Michael Efroimsky and Peter Goldreich (2004) "Gauge symmetry of the N-body problem of Celestial Mechanics." Astronomy and Astrophysics, Vol. 415, pp. 1187–1199. (2004)
- CODDINGTON, Earl.; LEVINSON, Norman. Theory of Ordinary Differential Equations. New York: McGraw-Hill, 1955. Dostupné online.
- BOYCE, W. E.; DIPRIMA, R. C. Elementary Differential Equations and Boundary Value Problems 8th Edition. [s.l.]: Wiley Interscience, 1965., stránky 186-192, 237-241
- TESCHL, Gerald. Ordinary Differential Equations and Dynamical Systems. Providence: American Mathematical Society Dostupné online.[nedostupný zdroj]
Externí odkazy
[editovat | editovat zdroj]- Paul Dawkins Lamar University: Online poznámky a důkaz.
- Stránka PlanetMath Archivováno 12. 12. 2009 na Wayback Machine..
- Motivace metody nebeskou mechanikou, (původní URL)